Incommensurate dynamics of resonant breathers in Josephson junction ladders 
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We present theoretical and experimental studies of resonant localized resistive states in a Josephson 
junction ladder. These complex breather states are obtained by tuning the breather frequency into 
the upper band of linear electromagnetic oscillations of the ladder. Their prominent feature is the ap- 
pearance of resonant steps in the current-voltage {I-V) characteristics. We have found the resonant 
breather-like states displaying incommensurate dynamics. Numerical simulations show that these 
incommensurate resonant breathers persist for very low values of damping. Qualitatively similar 
incommensurate breather states are observed in experiments performed with Nb-based Josephson 
ladders. We explain the appearance of these states with the help of resonance-induced hysteresis 
features in the I-V dependence. 
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A lot of attention has been recently devoted to the 
theoretical and experimental studu-qf intrinsic dynamic 
localized statesr-iu various physicaltru, chemical and bio- 
logical systemstra. These peculiar inhomogeneous states, 
called discrete breathers, appear in spatially homoge- 
neous (Jjlamiltonian or dissipative) nonlinear discrete 
systemsEl€l. It was found that these states manifest 
themselves in the spectra of molecules and solids, and 
can determine the energy transfer in complex biolog- 
ical systems!! Discrete breathers have been observed 
e.g. in weakly coupled optical waveguides, in anti- 
ferromagnetitspl^ds and in systems of interacting Joseph- 
son junctions! 
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FIG. 1. Josephson junction ladder. Crosses mark the in- 
dividual junctions. Arrows indicate the direction of external 
current flow (dc bias 7). The black circles indicate the junc- 
tions that are in the whirling (resistive) state. The breather 
state displaying up-down symmetry is shown, 

Especially in the latter example, namely homoge- 
neously dc driven Josephson junction ladders (J JLs) (see 
Fig. 1), an enormous diversity of breather states has been 
observed. The JJL consists of small Josephson junctions 
of two types, "horizontal" and "vertical" ones, which are, 
respectively arranged perpendicular and parallel to the 
applied dc bias current 7. Breather states in a JJL are 
characterized by spatially localized voltage patterns as 



a few Josephson junctions are in the resistive (whirling) 
states (the black circles in Fig. 1) while the rest of the 
JJL junctions are in the superconducting state. The dy- 
namics of breather states crucially depends on two pa- 
rameters: the anisotropy r\ — j^ 2 -, where I c h and J c y 
are respectively the critical currents of horizontal and 
vertical junctions, and the discreteness parameter (nor- 
malized inductance of the cell), [3l- 

The complete dynamics of JJLs is determined by the 
time dependent Josephson phases 0| (t) (vertical), 
(lower horizontal), and <${t) (upper horizontal), and is 
governed by the Kirchhoff's laws and the flux quantiza- 
tion law. By making use of the resistively shunted model 
for Josephson iqnctions the following set of equations has 
been obtainedli3: 

4>l + <*il + sinC = 7 + -^-(AC + V^_ x - V#U) 

PL 

4> h n + afa + sin <fa = --Uv< + <ft ~ 4>n) (1) 

VPL 

It + at + sin^ = -Uv^ + <t> h n - fa) , 

where we use the notations A/ n = f n -i — 2fn + fn+i and 
V/ n = fn+i — fn- Here, the unit of time is the inverse of 
the plasma frequency lo Pi and a is the damping constant 
for a single junction. The dc voltage drop across a single 
junction is defined through a time average < 4>{t) >. Of 
crucial importance for the rich variety of breather states 
is the presence of a nonzero dissipation in JJLs. This 
dissipation ensures that breathers can survive resonant 
interaction with extended linear electromagnetic waves 
(EWs), while their dynamic complexity may actually in- 
crease at the same time. Indeed, these complex resonant 
breather states have been predicted in [Ref. ^ , experi- 
mentally observed in [Ref. EEL and analyzed in detail in 
[Ref. [l4| . The condition of the appearance of such a res- 
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onant breather state is the matching of the breather fre- 
quency or its higher harmonics with the frequency Lu(q) 
of one of the cavity modes of the ladder, co(q) = mO. 
The breather frequency f2 is determined by the smallest 
nonzero dc voltage drop in the JJL, and the spectrum of 
cavity modes oj{q) has been derived previously using the 
linearization, nf. the Eqs. ([!]) around the superconducting 
ground stateoEfl: 

«4 = F ± VF 2 -G , 

F = \ + ^- + - 7 2 + ^-(1 - cosg) , (2) 
2 PlV 2 fi L 

G = (1 + -^-)\fl^f + ^-(1 - cosg) . 

For a finite size ladder with open boundary conditions, 
the spectrum of linear waves is discrete and characterized 
by the following choice of allowed wave number values: 

1l = J^r[ 1 = 0,1,2...,^, (3) 

where N is the number of cells. Note here, that the res- 
onant breather state is still a strictly periodic solution 
in time. Its resonant character is revealed by large am- 
plitude librations of the Josephson phases in the whole 
JJL. Resonant breather states are manifested by reso- 
nant steps in the current-voltage (I-V) characteristics, 
i.e. by a region where the voltage saturates at a fixed 
value while the bias 7 is increased. 

The presence of the resonant breather state leads to 
the appearance of hysteresis in the I-V curve. Conse- 
quently there exist particular ranges of 7 where both res- 
onant and non-resonant breathers can be stable. How- 
ever, these two types of breathers are characterized by 
different (incommensurate) frequencies fi, and more- 
over, the frequencies depend in a different manner on 
7. Thus, the question naturally appears: Is it possi- 
ble for a breather state to be both resonant and non- 
resonant at the same time in a JJL? This is important be- 
cause nondissipative systems support only time-periodic 
breathers, while dissipative systems typically need an ex- 
ternal time-dependent drive to generate quasiperiodic dy- 
namics. The idea here is that the spatial nonuniformity 
of the states (localized breathers) acts in a similar way 
to external time-dependent drives, so that even a dc bias 
is enough to obtain quasiperiodic (in time) states. 

Here we will consider a region of parameters where 
the breather strongly interacts with the EWs. We find 
that for low values of damping a a novel breather-like 
state with incommensurate (aperiodic) dynamics per- 
sists. This incommensurate resonant state can be viewed 
as the mixture of resonant and non-resonant breathers 
in the JJL. Moreover, we also report on the experimen- 
tal observation of this peculiar breather state in an ex- 
tremely underdamped JJL. The incommensurate reso- 
nant breather-like state leads to two successive voltage 
jumps in the I-V curve. 



To study the breather dynamics we perform direct nu- 
merical simulations of the set of equations ([[]). To es- 
tablish a large dc current region where a strong resonant 
interaction of the breather state with the EWs occurs 
we choose the discreteness parameter (3l = 0.378. The 
ladder consists of N = 10 cells and the anisotropy param- 
eter 77 = 0.49. We impose open boundary conditions 
and use 4th order Runge-Kutta method for integration. 
The initial value of the dc bias was 7 = 0.8. We select 
proper initial conditions that lead to the relaxation of 
the system into a particular breather state of up-down 
symmetry with one resistive vertical junction, as shown 
in Fig. 1. Similar to our previous workEJ, after a waiting 
time of 500 time units we use the next 500 time units to 
calculate the time averaged characteristics of the state. 
Then the dc bias 7 was decreased (increased) by a step of 
A7 = 0.001 and the procedure was repeated. We checked 
that our results do not change upon further increase of 
the waiting time. The I — V characteristics, that is the 
dependence of the averaged voltage drop across the resis- 
tive Josephson junction on the dc bias, was monitored. 
The simulations were carried out for two different values 
of damping a — 0.1 (see Fig. 2) and a — 0.025 (see 
Fig. 3). 
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FIG. 2. I - V characteristics for a = 0.1, (3 L = 0.378, 
77 = 0.49. Arrows show the dc bias current increase (decrease). 
The dc voltage pattern of the breather state is presented by 
the black circles. 



In the case of a moderate value of damping a = 0. 1 we 
obtain a non-resonant breather state for dc bias current 
values down to 7 > 0.72. This breather state displays 
a complete up-down symmetry, $(t) = — (t ) . The 
I-V curve is linear in this regime. In the region of mod- 
erate dc current bias values resonant steps appear in the 
I-V curve (Fig. 2). The presence of these steps indi- 
cates a strong resonant interaction between the breather 
state and EWs as the condition w+(<7z) = 2£! is satis- 
fied. The observed steps correspond to resonances with 
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cavity modes with particular values of I = 2,4,6,8. 
As the JJL is biased at the steps, a resonant type of 
the breather state is realized. This state conserves the 
structure of the dc voltage pattern but is characterized 
by a large ac voltage component of the superconducting 
Josephson junctions. Independent of the observed dy- 
namical complexity, all breather states (resonant or non- 
resonant type) found for this value of damping a display 
a time-periodic behavior. 

However, the situation changes drastically as we lower 
the dissipation down to a = 0.025. First, the 
breather state strongly interacts with EWs if the condi- 
tion u>(qi) = is matched. Secondly, we observe only a 
single resonant step corresponding to the particular cav- 
ity mode (1 = 2) excitation. Under these conditions we 
founcLihat the switching from the non-resonant breather 
stateE3 appearing for large values of dc voltages, to the 
resonant breather occurs in the form of two successive 
voltage jumps (see, Fig. 3a). 
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novel incommensurate breather-like state can be viewed 
as a coexistence of resonant and non-resonant junction 
behavior in the JJL. We note here that similar aperiodic 
breather states were found in a wide range of parameters 
P L ~ 1 and r\ < 1. 

By making use of a dc analysisalij, i.e. neglecting ac 
components of the Josephson current, we find that the 
voltage jumps AVh and AVh of the non-resonant and 
resonant horizontal junctions have to satisfy the relation- 
ship: 



AV h = 



AV h 
I + 277 



(4) 



Moreover, the presence of the incommensurate resonant 
breather state is characterized by a lower resistance oc 
(a(2r/+ 1)) _1 as compared to the non-resonant case. In- 
deed, our simulations show that the dc voltage jump of 
the lower horizontal junctions is approximately two times 
smaller than the corresponding dc voltage jump of the 
upper horizontal junctions, in full accordance with the 
dc analysis. 

In order to characterize an incommensurate breather 
state more precisely we calculate the Fourier transform of 
the time dependent voltage 4>\(t) of the vertical junction 
(see Fig. 4), which is in the superconducting state and 
the closest one to the core of the breather. 
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FIG. 3. I -V characteristics for a = 0.025, 0l = 0.378, 
f] = 0.49. Arrows indicate the voltage jumps and show the 
decrease of the dc bias current, a) Result for voltage drop 
across the vertical junction. The size of the black circles in- 
dicates the relative values of dc voltage drops, b) Result for 
twice the voltage drop across the upper (solid line) and lower 
(dashed line) horizontal junctions. 



Thus, in the particular region of the dc bias current 
(0.31 < 7 < 0.38) the up-down symmetry of the dc 
voltage pattern is violated and a novel resonant breather- 
like state with incommensurate (aperiodic) dynamics ap- 
pears. Indeed, as the system jumps to this state there 
are at least two incommensurate frequencies, namely 
ill — 3.4 and O2 ~ 5.3. Moreover, with decreas- 
ing dc bias the value of fii, practically does not change, 
in contrast to O2. A closer look at the I-V curves of 
the horizontal junctions (Fig. 3b) also shows that the 
upper (lower) horizontal junctions display the resonant 
(non-resonant) behavior. Thus we can conclude that this 
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FIG. 4. The Fourier transform (FT) of the ac voltage <j>X(t). 
Three different values of dc bias current are shown and the 
logarithmic scale for FT was used. 



If the breather state is the periodic one , i. e. for large 
(Fig. 4a) and small (Fig. 4c) values of dc bias current, 
the Fourier transform contains the breather frequency f2 
and higher harmonics. However, in the intermediate re- 
gion of the dc bias as the breather displays incommensu- 
rate dynamics, the Fourier transform becomes more rich 
(Fig. 4b). We observe two different harmonics fii and 
O2, and the various combination frequencies mQi + n£l2, 
where n and m are integers. 
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The incommensurate breather-like states analyzed 
above were initially observed in experiments which we 
performed using Nb-based Josephson junction ladders. 
The experimental layout is similar to the one presented 
in [Ref. [| and corresponds to the schematic shown in 
Fig. 1. The measurements were performed in liquid he- 
lium at a temperature of T = 4.2 K. The data we present 
were taken from a 10 cell open-ended Josephson ladder of 
anisotropy r\ = 0.49. By chosing a relatively low critical 
current density of j c m 100 A/cm 2 for the samples we 
achieved a discreteness parameter of (3l — 0.62, which 
makes it possible for a breather to be strongly resonant 
with EWs. At the same time, the low damping value of 
a = 0.025 was sustained. 

The breather state was created artificially by initially 
applying a local bias current 7/ to one of the vertical 
junctions, and thus, forcing the switching of this junction 
into the resistive state. After that we decreased 7; and 
simultaneously increased the homogeneous bias current 
7. This procedure, described in detail in [Ref. 0], finally 
yields the desired breather state in the presence of an 
uniformly supplied bias current. In this paper we focused 
on an up-down symmetric breather state with dc voltage 
pattern shown in Fig. 1. The I-V curves were obtained by 
sweeping the bias current 7 and measuring the respective 
voltage drops across the junctions. For convenience, we 
plot normalized voltage v — V/Vq with Vq = = 
84.3 [j,V, and normalized bias current 7 = I ex t/IcV with 
I c = 20.5 fiA. 

For the sample #1 the breather state was created at a 
dc bias current 7 = 0.38. The acquired I-V curve, as 
shown in Fig. 5, displays a linear, non-resonant behav- 
ior until 7 = 0.35 (see Fig. 5). The breather state 
is symmetric, and the absolute values of the horizon- 
tal junction voltages \Vh\ and \Vh\ are just one half of 
the vertical junction voltage, \V V \. As the dc bias is de- 
creased further, the switching to the resonant breather 
state takes place. This switching occurs through an in- 
termediate state (see Fig. 5a). Indeed, in this region of 
0.32 < 7 < 0.35 we observe the symmetry-broken state 
with \ Vh\ 7^ \ Vh\- Moreover, the ratio of voltages |Vh|/|14| 
continuously changes as 7 decreases (see Fig. 5b) . For the 
values of 7 < 0.32 the symmetric breather state recovers 
but is now a resonant breather state. 

For the sample #2 we observe a similar intermediate 
state appearing in the process of the switching from the 
non-resonant to the resonant breather state (see Fig. 6). 
In this experiment carried out in the presence of a small 
magnetic field (the normalized magnetic flux per cell 
was less than 0.1), the I-V curve of the non-resonant 
(lower) horizontal junction (see Fig. 6b) displays a volt- 
age jump similar to the one found in the numerical stud- 
ies (see Fig. 3). The ratio of the voltages and the cor- 
responding rotation frequencies of the horizontal (upper 
and lower) junctions increases as 7 decreases. The mea- 
sured I-V curves in both cases qualitatively look similar 
to the current-voltage characteristics obtained by numer- 
ical simulations of (1). Thus, the predicted breather state 



displaying aperiodic dynamics is clearly observed (com- 
pare Figs. 3 and 5,6). 
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FIG. 5. Measured I-V curves (sample #1) of an initially 
symmetric breather state created at a dc bias of 7 = 0.38: (a) 
voltage drop across vertical junction, (b) voltage drops across 
upper (thin solid line) and lower (dashed line) left horizontal 
junctions. The breather state with incommensurate dynamics 
is formed for 0.32 < 7 < 0.35. The fine structure superim- 
posed on the curve is due to reduced experimental resolution 
in these measurements. 



However, there are several differences between the ex- 
perimental observations and the numerical analysis of 
(1). Firstly, in contrast to numerics we did not observe 
the voltage jump of the non-resonant horizontal junction 
in the measurement shown in Fig. 5, and secondly, in 
Fig. 6 we observe just a tiny voltage difference between 
upper and lower horizontal junctions in the bias region 
above the incommensurate state. This particular behav- 
ior is not yet explained. We want to stress that the in- 
commensurate breather state was observed under various 
experimental conditions, which were partly influenced by 
magnetic fields in the vicinity of the studied sample. 

In conclusion, we found an inhomogeneous resistive 
state showing aperiodic dynamical behavior by making 
use of a numerical analysis of (1). Similar incommen- 
surate breather states were observed experimentally in 
niobium-based underdamped Josephson junction ladders. 
The incommensurate breather states are characterized by 
a resonant behavior of some horizontal junctions and a 
non-resonant behavior of the other horizontal junctions. 
Such states appear in the process of a switching from 
symmetric non-resonant to symmetric resonant breather 
state. The origin of these particular breather states is a 
strong resonant excitation of EWs by the breather state 
leading to the resonance-induced hysteresis in the I-V 
curve. Thus, we argue that the aperiodic inhomogeneous 



4 



dynamic states can present a generic feature of dissipa- 
tive (but extremely underdamped) lattices. 
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FIG. 6. Measured I-V curves (sample #2) of an initially 
symmetric breather state created at a dc bias of 7 = 0.4: (a) 
voltage drop across vertical junction, (b) voltage drops across 
upper (thin solid line) and lower (dashed line) left horizontal 
junctions. The breather state with incommensurate voltages 
is formed for 0.34 < 7 < 0.37. 
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